NONTHERMAL HEATING 

: Wv 7 ^ 77 : IN 'THE -7 7 ‘ ^ 

TWO-FLUID SOLAR-WIND MODEL 








R. E. HARTLE ' 
AARON BARNES 


7 '.. ' !| 




\ \ 


APRIL 1970 


V ! " 

J,-l. ■ 


) : ,7i * 

V 7 , 

?V ^ 

■ ~ > 


GSFC 


GODDARD SPACE FLIGHT CENTER 

I . GREENBELT, MARYLAND 


, / 

■ ' ‘ 1 


NONTHERMAL HEATING IN THE TWO -FLU ID 
SOLAR-WIND MODEL 

by 

R. E. Hartle 

Laboratory for Planetary Atmospheres 
NASA, Goddard Space Flight Center 
Greenbelt, Maryland 20771 

Aaron Barnes 
Space Sciences Division 
NASA, Ames Research Center 
Moffett Field, California 94035 

ABSTRACT 

The two-fluid solar-wind model predicted values for wind 
speed and proton temperature that are lower than average ob- 
served values; however, the predictions are consistent with 
the empirical relation of Burlaga and Ogilvie in which the 
square root of the proton temperature is proportional to the 
wind speed. Since observed values of wind speed and proton 
temperature are typically higher than those of the two-fluid 
model, it is concluded that energy must be supplied by dissipa- 
tion of nonthermal energy from an external source. The two- 
fluid model is therefore extended by including a<i hoc an energy 
source in the proton heat equation whose strength and spatial 
distribution is varied to determine the general requirements 
of such a source. The main features of our results are: 

Heat deposited inside the heliocentric radius r - 4 Rq results 
in a significant increase in wind speed with negligible in- 
crease in proton temperature. Depositing heat over the range 
r ~ 25 Rq results in a large increase in proton temperature and 
negligible increase in wind speed . By depositing heat over the 
extended range 2 R© ~ r ~ 25 R©, solar -wind speeds and proton 
temperatures can be brought in direct correspondence with the 
empirical results of Burlaga and Ogilvie . Based on this model 
we conclude that primary energy deposition should take place 
inside r - 25 R©. 



1. INTRODUCTION 


In this paper we describe a simple two-fluid model of solar- 
wind flow which includes dissipation of non thermal energy as well 
as thermal conduction. From its beginning, the theory of the 
solar wind has been concerned with the problem of energy trans- 
port in the solar corona [see reviews by Parker , 1969 , 1965a , 

1963] . The coronal gas is fairly tightly bound in the sun’s 
gravitational field , so that if there is a solar wind , there 
must exist a heat source strong enough to lift the gas out of 
solar gravity. This heating must extend out from the sun at 
least as far as the critical distance (the distance at which the 
solar-wind flow becomes supersonic) . 

The mechanisms which provide the necessary heating are not 
yet understood in quantitative detail. Parker [1963, 1965a, b] 
has discussed in detail the possibility that the heat transport 
of the wind is entirely due to thermal conduction. He found 
that purely conductive models imply lower wind velocity and 
higher density than are observed at the orbit of earth (in fact, 

g o 

f or a base temperature of 2 x 10 K, the observed coronal density 
9 -3 

of ~10 cm implies subsonic expansion) , a conclusion supported 
by the numerical calculations of Noble and Scarf [1963] and 
Whang and Chang Cl965] . Parker suggested that at least two 
other effects may be important in solar -wind energetics : non- 
conduct ive heating , probably due to wave dissipation , and cutoff 
of thermal conduction, due either to spiraling of the interplanetary 
magnetic field or to some collisionless mechanism. Extended heating 
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in one-fluid models has been discussed in the polytrope ap- 
proximation by Parker [1963], Analytical discussions of 
one-fluid heating from several viewpoints have been given by 
Konyukov [l967a, b, 1968]. Konyukov' s arguments are not 
completely convincing, however, especially as he incorrectly 
concludes that a purely conductive one-fluid model is unstable 
[Konyukov, 1968], 

A more refined solar-wind model, where the only heat 
transport is classical heat conduction, is the two-fluid model 
[ Sturrock and Hartle , 1966, Hartle and Sturrock , 1968], in 
which the distinct transport properties of electrons and pro- 
tons are considered. This model predicts a flow speed of 

3 o 

~250 km/sec and a proton temperature of ~4 x lO' 3 °K at the 
orbit of earth, considerably smaller than the typically ob- 
served values of ~ 350 km/sec and ~ 5:. ±lf0^ °K [ Neugebauer and 
Snyder , Wolfe et al , 196J6; Axford, 1968, Hundhausen , 1968]; 

in addition, the two-fluid model predicts an electron temperature 
5 o 

of ~3. 5 x 10 K at 1 AU , larger than the observed value of 
5 o 

~2 x 10 K. It has been proposed that these deficiencies might 
be remedied by modifying this purely conductive model to include 
the possibility of heating by the dissipation of some nonthermal 
energy supply; several possible mechanisms have been discussed 
in the literature [ Coleman , 1968; Barnes , 1969; Fredricks , 1969; 
Jokipii and Davis , 1969]. Recently, Hundhausen [1969] has sug- 
gested that inhibition of thermal conduction, rather than dis- 
sipation of nonthermal energy from some external source, may 
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account for the discrepancies between quiet-day observations 
and the predictions of the two-fluid model . Such inhibition 
of thermal conduction might be either a purely geometrical 
effect, associated with the spiraling of the interplanetary 
magnetic field [ Brandt , et al, 1969; Gentry and Hundhausen , 

1969; Urch , 1969], or a collisionless plasma effect, perhaps 
of the sort discussed by Fors lund [1970]. 

Nevertheless, the observed solar wind often blows much 
faster than the quiet-day wind for periods of several days 
(e.g. , see Neugebauer and Snyder , 19®(j3&;) , so that the wind 
convects much more energy at some times than at others. This 
variability of the energy flux is almost certainly associated 
with the dissipation of varying amounts of nonthermal energy 
[ Parker , 1969, 1965a], This viewpoint is consistent with the 
strong observational correlation between the flow speed v^.. and 
proton temperature T of the wind at 1 AU [Burlaga and Ogilvie, 

PE 

1970; Robbins et al,1970]. In particular, Burlaga and Ogilvie find 
the empirical relation T^ 2 = (.036 ± .003) v £ - (5.54 * 1.50), 
where the proton temperature T _ is in units of 10^ °K and the 

PiL 

flow speed v E is in km/sec. (subscript E refers to Earth's orbit 
throughout) . This remarkable relation is valid for a wide variety 
of solar-wind conditions 1 , ranging from quiet to highly disturbed 
times. There is also evidence for a weaker correlation between 
density and temperature [ Robbins et al , 1970] , especially for temp- 
eratures lower than 5 x 10 4 °K [ Hundhausen , private communication] . 

4 This relation is based on Explorer 34 plasma measurements made 
between June and December of 1967 , but seems to agree with other 
observations made at other times in the solar cycle . 
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However, the T pE ~v E correlation is much stronger [ Robbins et al , 
1970], and probably a more fundamental characteristic of solar- 
wind flow at the orbit of earth. 

Since disturbed wind conditions often persist for several 
days (a time comparable to the time of flow between the sun 
and the orbit of earth) , it seems likely that the Burlaga- 
Ogilvie T_„-v„ relation defines a continuum of average macroscopic 

pill ill — 

states of the solar wind,^ rather than mere correlations in the 
microstructure of the wind; presumably, high velocities (and 
temperatures) indicate a state of the wind that has undergone 
extended heating. If this interpretation of the T pE -Vg relation 
is correct, it has important implications for theoretical models 
of the wind. For example, most one-fluid models [ Noble and Scarf , 
1963; Whang and Chang , 1965; Gentry and Hundhausen , 1969, Urch, 
1969] give values of flow speed and temperature that separately 
lie in the observed range, but do not satisfy the T pE “ v E relation; 
accordingly, we conclude that none of these models describes an 
observed average state of the wind. On the other hand, the two- 
fluid model is consistent with extrapolation of the T p E -v E rela- 
tion to low temperature and flow speed. 

If the Burlaga-Ogilvie T -v relation is a fundamental 

pE E 

property of macroscopic solar-wind flow, it must be explainable 
by a sufficiently realistic model of the wind. Obviously, such 
a model must not simply predict a flow speed and proton tempera- 
ture at (say) a heliocentric distance of 1 AU but rather it should 
predict a continuous range of T and v consistent with the 
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empirical T p E -v E relation. In the present paper, we discuss 

a simple model which is consistent with the empirical T pE ~v E 

relation ; although this model does not actually predict the 

T „ — v_ relation, we believe that it indicates some of the 
piL h 

features that a predictive model must have. Our present model 

is similar to the two-fluid model in the form presented by 

Hartle and Sturrock [1968], but includes heating of the ions 

(but not electrons) by dissipation of nonthermal energy from 

2 

an external source. We do not specify the mechanism that 
supplies this nonthermal energy, but assume ad hoc an energy 
dissipation function which permits us to vary the strength and 
spatial distribution of the heating in a general way. 

We solve the dynamical equations of solar-wind flow for 
various shapes and strengths of the energy deposition function. 

We find that heat deposited in the region of subsonic flow 
mainly raises the flow speed v E at 1 AU, and heat deposited in 
the region of supersonic flow mainly raises the proton tempera- 

3 

ture T _ at 1 AU. By discarding solutions not consistent with 

p£j 

the Burlaga-Ogilvie T pE ~v E relation, we obtaiin a series of solu- 
tions which indicate that the external heat source has the following 
property: the primary heating occurs over an extended region, 
possibly as far as 25 R 0 from the sun; also, weaker heating may 


2 

This model is fluid theory, as opposed to exospheric theories 
(e.g. , Hollweg [1970] , Eviatar and Schulz [1970 ] , Jockers [j.969]) . 


3 

This effect is not surprising , since the addition of heat to 
subsonic flow increases Mach number, while addition of heat to 
supersonic flow reduces Mach number [Shapiro, 1953] . 
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occur much farther out (~ 1 AU) , but has only a secondary 
influence, affecting the temperature but not the flow speed. 

In the limit of zero energy deposition our solutions reduce 

4 

to the basic two-fluid model. Although this limit may not be 
an attainable state of solar wind flow, because the large ratio 
of electron to proton temperature may cause instability [ Fredricks , 
1969, Forslund , 1970], it is reasonable to use the two-fluid' 
equations to investigate the effects of nonthermal energy deposi- 
tion for two reasons: first, the addition of modest amounts of 
heat raises the proton temperature high enough that the strongest 
instabilities considered by Fredricks and Forslund do not occur; 
second, as mentioned before, extrapolation of the empirical 
T „-v_ relation to low velocity and temperature is consistent 

Pill J Li 

with the basic two-fluid model, suggesting that it is a reason- 
able mathematical zero-heating limit of the models with finite 
heating. 

Hundhausen [1969] has emphasized that the two-fluid model 
implies a thermal conduction flux considerably higher than that 
deduced from electron measurement at 1 AU [ Montgomery et al . , 

1968]. He suggested that this may mean that the transport co- 
efficients of the two-fluid model become: invalid at some point 
in the flow; this possibility had already been mentioned briefly 
by Hartle and Sturrock [1968] and has been discussed in consider- 
able detail by Forslund [1970] . It is clear that inhibition of 


4 

One cannot produce the T _,-v relation by variation of the base 
parameters of the two-fluid model. This was shown by numerical 
computation, and also may be seen from the approximate equations 
<§.l) 1 0feH§g4le HrifloS^urrock [1968] . 
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thermal conduction, by whatever mechanism, is an effect which 
should be taken into account in a realistic model of the minimum 
solar wind. However, it is doubtful that inhibition of thermal 
conduction is related to the wide variation in energy flux 
associated with various macroscopic states of the solar wind. 

This kind of effect has been omitted from the present model, 
and probably does not modify our qualitative conclusions in a 
significant way. 

Hundhausen [1969] has also raised the point that the two- 
fluid model predicts a total energy flux at 1 AU which is larger 
than the observed quiet-day flux by a factor ~2. This may not 
be a fundamental difficulty of the two-fluid model, for the 
total efflux of energy in the quiet-day wind is a fairly small 
fraction (~ 10%) of the heat required to drive the wind (most 
of this heat goes into work of lifting the plasma out of solar 
gravity) , so that a small error in describing the heating process 
may cause a rather large error in the net energy flux in the 
distant flow. Further, the observational correlation between 
temperature and flow speed is much stronger than that between 
either of these quantities and density, so that density (and 
hence energy flux, whose convective part is proportional to 
density) may not be a critical test of a solar-wind model £cf. 
also Parker, 1965a]. 
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2. BASIC EQUATIONS AND THEIR SOLUTION 


To make the calculations tractable we assume that the 
solar wind consists of protons and electrons and undergoes 
steady, radial, spherically symmetric expansion characterized 
by a proton number density n and flow speed v (equal to the 
electron density and flow speed) , proton temperature T and 

ir 

electron temperature T . Effects of solar rotation, magnetic 
field, temperature anisotropy and viscosity are ignored. The 
equations that govern the flow are the continuity equation 


n v r = J = constant 


(1) 


(r is the heliocentric distance) , the momentum equation 

dv d 


n m p V dr 


G m n 

ar [nk(T e + V 3 7T- 


( 2 ) 


and the electron and proton energy equations 


3 . , _ dn 

J nvk 3T- vkT edF 


1* /r 2 K 
r 2 dr ^ e dr J 


S' E 


v_ nk(T - T ) + « 


(3) 


dT 


j? n v k r 
2 dr 


p - v k T as. 

p dr 


_ L_ d ( r 2 K dT P "\ = 

r 2 dr \ r p dr J 


3 

2 


v E n k(T e 


V 


+ P, 


(the P’s are unspecified collisionless heating terms) . 


(4) 
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Here nip and m e are, respectively, the proton and electron 
masses, G is the gravitational constant, M 0 is the mass of 
the sun, k is Boltzmann’s constant, 


n k 2 T 


1C = 3.I6 

e m v 

e e 


( 5 ) 


and 


tc 

P 


3.9 


n k 2 T 

m v 
P P 


( 6 ) 


are the electron and proton collisional thermal conductivities, 


2m 

v = — — v 

E m p e 


(7) 


is the collisional electron-proton energy exchange rate, 


and 


_ 4(2rr) 1//2 n e 4 \ 

v e Trr—r^ri 


3 m. 


(k V 


v = 

P 


A 1/2 4 

4tt n e \ 


( 8 ) 


(9) 
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are the electron and proton collision frequencies, and 



( 10 ) 


is the Coulomb logarithm. The expressions for the collisional 
transport terms are taken from the review article by Braginskii 
[1965 '0. 

Equations (1-10) are the same as the basic equations used 
by Hartle and Sturrock [1968; hereafter called HS>] except for 
the collisionless heating terms P g and P^ in equations (3-4) . 

It is proper to treat the collisionless heating by additive 
terms, as corrections to the particle kinetic equations due to 
collisionless effects (wav e- part icle interactions) are additive 
(e.g., see Rogister and Oberman , [1968] ) . In principle the 
momentum equation (2) (but, obviously, not the continuity 
equation (l))has an additional term due to collisionless effects, 
a pressure-like term analogous to radiation pressure. How- 
ever, one may argue on physical grounds that this collisionless 
acceleration is not a major effect, and that the collisionless 
processes influence acceleration only through modifying proton 
and electron temperatures according to the heat equation. For 


suppose that the collisionless heating is due to some kind of 
dissipation of fluctuations of energy density ^ \ which move 
with speed u relative to the wind, and that u is comparable to 


the proton thermal speed or Alfvdn speed; 


of the fluctuations is of order 



the momentum density 
The corresponding 


gas energy density appearing in the heat equations is of the 
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order of the pressure P = nk(T + T ) , and the related gas 

e p 

momentum density is of order P/v + m nv . We estimate the 

P 

importance of the collisionless terms in the heat equat ions 
by the ratio n of the energy flux of the fluctuations to the 
gas energy flux appearing in the heat equations , 

T| ~ ^ U ' Pv~ V ' ^' = p (1 + U / V ) * 


If n is of order 0.1 or larger, we expect collisionless heating 
to be important in the dynamics of the wind. Similarly, we 
estimate the importance of collisionless terms in the momentum 
equation by the ratio 


u 




u+v 


(u+v ) 2 * 2 — 

P+m nv (1+u/v) (1+m nv /P) 
P P 


It is clear that n « n both in the subsonic-flow region, where 
2 

u » v and m^nv « P, and in the supersonic-flow region, where 

2 

u « v and m nv » P. Hence collisionless processes are much 
P 

less significant for acceleration than for heating (except 
possibly in the small region of subsonic-supersonic transition) , 
and terms describing collisionless acceleration need not appear 
in equation (2) . 

The solutions of the dynamical equations (1-4) depend on 

the forms of the collisionless heating terms jp and p . In the 

e p 

present work we assume that P e = 0 and 


p p " D o <£“> ex P 
* o 


(r/R - a)' 


~T 


( 11 ) 


11 



D q is a measure of the strength of collisionless heating 
(D q = 0 corresponds to the bare two-fluid equations of HS) , 
a is the heliocentric distance (in R 0 ) of maximum energy 
deposition, b is the spread of the energy deposition, and 
n Q is the proton number density at the base of the model 
r = r . The total energy deposition rate by collisionless 
processes in the wind between the base of the corona and 
heliocentric distance r is 


AE(r) 


4nD o i* exp 

R 0 ° 


(p/R Q -a) 2 

„ 

b 2 



( 12 ) 


which must be computed numerically from the solutions of the 
dynamical equations. We emphasize that equations (11) and (12) 
are not realistic descriptions of the energy deposition; the 
form (11) was chosen because the exponential term permits 
straightforward investigation of the qualitative effects of heat 
sources of various strengths, locations and spreads, and the 
factor (n/n Q ) permits a simple dependence of the heating on the 
amount of matter present . 

Except for the specification of the heating strength D q , 
its position a, and its spread b, the boundary conditions are 
the same as in HS. We specify the base density and temperatures 
n o , T eQ , T p 0 > require that T 0 and T p tend to zero as r goes to 
infinity, and that the solar -wind velocity should pass continuously 
through a subsonic-supersonic transition. 

The dynamical equations (1-4) are solved numerically by the 
iterative technique of HS. We begin by guessing the temperature 
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profiles T (r) and T (r) , and use them to solve the momentum 
e p 

equat ion (2) (from which the density has been eliminated via 
(1) ) . There is a unique solution of equation (2) which exhibits 
the required subsonic-supersonic transition; this velocity pro- 
file v (r) , together with the base density n Q , give the density 
profile n (r) according to equat ion (1). 

We then use the values of v (r) so computed , together with 
the assumed T (r) , to solve the electron heat equation (3) for 

r 

T e (r), subject to the boundary conditions on T g at r = and 
r = ce. The outer boundary condition is handled by requiring 

o /7 

that T_(r) exhibit the asymptotic dependence T (r) « r 

at large r; this asymptotic form obtains because the conduction 

term dominates the electron heat equation at large r. Having 

calculated n(r), v(r) and T (r) , we use these data to solve 

the proton heat equation (4) for T (r) , subject to the boundary 

P 

condition on T p at r = t Q . This equation is essentially a first- 
order equation beyond a few solar radii, because the heat exchange 
term dominates the proton conduction term. In this case, 

— /7 

T (r) cc r" 7 as r -* co because of the dominance of the energy 
P 

—2 /7 

exchange term and the r dependence of T . 

This sequence of calculations, which began with assumed 

values of T (r) and T (r) , has yielded new functions T and T 
0 P © P 

satisfying the boundary conditions . These may be used to solve 
the equation of motion again and so repeat the previous cycle . 

This cyclic process is repeated until the changes in the variables 
from one cycle to the next are so small as to be judged in- 
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significant. To this approximation, all equations and boundary 
conditions are satisfied. Further details of the solution pro- 
cedure are discussed in HS . 

The introduction of the heating term (11) into the two-fluid 
model calculation introduces essentially no new element into the 
computation procedure. However, in order to get good first 
"guesses" of the temperature profiles, it is convenient to generate 
a series of models, beginning with a bare two-fluid model 
(D q = 0), then gradually raising D q in steps, using in each 
step the temperature profiles computed in the previous step 
as the basis for choosing the first trail temperature profiles 
for the iterative solution of the dynamical equations. 
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3. ANALYTICAL RESULTS 


The principle model of HS, chosen on the basis of a 

close fit to the observed coronal density profile of Blackwell 

[1956] , required n = 3 x lO^cm -3 and T = T = 2 x 10 6 °K 
’ ^ o eo po 

at a base radius r = R_ . We have used this model as a point 
of departure to determine the effects of collisionless proton 
heating. However, we have taken a base radius r Q of 2R Q to 
avoid the necessity of including the complications due to 
strong small scale magnetic fields and attendant non-radial 
expansion expected to be important inside this point. Then, 
for this case, the derived base values of HS we have used are: 
n Q = 1.5 x 10 6 cm" 3 , T 0Q = 1.5 x 10 6 °K, and T pQ - 1.2 x 10 6 °K. 

We have divided the analysis into two classes based on the 
parameters of the assumed density weighted gaussian heat source 
given by equation (11). In the first class, the half width b 
is held fixed while the central position of the gaussian 
a (a > 5) and the amplitude factor D Q are allowed to vary. In 
the second class, a is fixed at 2 while b and D Q are varied. 

The first class of models corresponds to energy deposition 
almost all of which takes place well above the base while the 
second class corresponds to heating over a finite range extending 
from the base. To determine the significance of each model in a 
given class, with heat source characterized by (a,b,D ), we compare 
the corresponding proton temperature T pE and wind speed v^, at 
Earth's orbit with the empirical relation of Burlaga and Ogilvie 
discussed above. 
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The analytical results of the first class of models are 

shown parametrically by the solid lines (and their extensions 

shown by dotted lines) in Figure 1 where the square root of 

the proton temperature at 1 AU in kilo-degrees Kelvin is 

plotted against the corresponding wind speed in km/sec. The 

resulting density n £ is indicated at selected points along the 

curve. The dashed lines are the empirical results of Burlaga 

1/2 

and Ogilvie, where the central dashed line, given by T p g 
0 . 036v„ - 5.54, represents the average values and the outer 

ill 

dashed lines, given by T = (0.036 + 0.003) v £ -(5.54 + 1.50), 

represent their variance. In this series, the half width b 

is fixed at 5 while each individual curve corresponds to a 

fixed central position a. Then, upon increasing the rate of 

energy deposition by increasing D Q from zero, for a given a 

and b, the solution points move upward from the basic solution 

of HS indicated at the apex of the curves. 

For completeness, the values of the amplitude factor D q 

corresponding to solution points designated by the density values 

and listed in an order starting with the basic HS solution are: 

for a = 6, D q x 10 8 = 1, 1.5; for a = 9, D q x 10 8 = 0.5, 1, 1.5, 

2.5; for a = 12, D q x 10 8 = 0.2, 0.6, 1, 1.5, 2; for a = 16 , 

D q x 10 8 = 0.3, 0.7, 1, 1.3; for a = 20 , D q x 10 8 = 0 .44 , 

0.97, 1.2; for a « 30, D q x 10 8 = 0.24, 0.58. For example, 

“8 

when a = 12 and v„ = 340 km/sec we note that D = 1,5 x 10 

JCj O 

The total energy deposition rate AE (°°) of equation (12) 
can be determined to a good approximation by 
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taking the difference between the streaming energy rate 

(at 1 AU) and the corresponding rate of the HS model without 

heating. 

Some of the results shown in Figure 1 are not consistent 

with the assumption of a local maxwellian velocity distribution 

implicitly required by the basic equations (2) through (4) . 

Local thermodynamic equilibrium is maintained when the proton 

momentum and energy relaxation rates v , given by equation 

P 

(9) , are greater than the characteristic expansion rate 
v = -(v/n) (dn/dr) . The analytical results lying in the 

cXp 

region where equilibrium is not expected to be maintained by 

Coulomb collisions, i.e. v < v , are indicated by dotted 

lines. Then, the results lying in the range of observation 

with wind speeds greater than about 320 km/sec and proton 

4 o 

temperatures greater than 3 x 10 K are not to be relied upon 

in the strictest sense. The results lying in the observed range 

with Vj, somewhat greater than 320 km/sec and T ^ somewhat greater 
4 o 

than 3 x 10 K are to be regarded as preliminary estimates 
to solutions of an appropriate kinetic equation required in this 
region. 

In every case shown in Figure 1, for sources of fixed spread 
b and position a, the solar wind speed v E and proton temperature 
T „ increase while the density n_ decreases when the heating rate 

pJa 

is increased (by increasing D Q ) . More importantly, this same 
variation is noted for those solution points consistent with 
the empirical relation of Burlaga and Ogilvie. In this case, 
the central position a is constrained to move inward toward the 
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sun while the heating rate is increased. That is, if D q is 
increased and a is decreased in a specified manner so that 
solution points lie on the central dashed line then v E and 
T „ increase while n„ decreases. In this case, we find that 
the amplitude D Q varies as the inverse square of the central 
position a over the range 275 km/sec < v^, < 400 km/sec. The 
decrease in density n £ is desirable since densities observed 
when v„ > 250 km/sec are usually lower than the 15 cm re- 

Jci 

suiting from the model of HS without heating. 

At this point we conclude that to attain solar wind 

speeds, say, greater than 300 km/sec while simultaneously 

maintaining proton temperatures low enough to agree with 

observation the primary energy deposition must take place 

within a heleocentric radius of about 25 R (see results for 

© 

a = 16, 12, and 9 in Figure 1). Heating outside this range 
produces primarily an increase in proton temperature T _ with 

pill 

negligible increase in wind speed v_ as indicated in Figure 1 

ill 

when a = 30. To further establish this conclusion, we have 

investigated the effects of varying the spread b, the base 

temperature Tp Q , and the base density n Q . 

Upon varying the half width b over the range 1 - b - 9, 

while keeping a and D 0 fixed, we find that solution points in 

the observational range move approximately along the corresponding 

curves of constant a in Figure 1. That is, increasing (decreasing) 

b corresponds to an increase (decrease) in v £ and Tp E as the 
1 /2 

point (Tp E , v £ ) moves along a curve of constant a. This 
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result follows since an increase or decrease in b implies an 
increase or decrease, respectively, in the rate of energy de- 


position. The corresponding variation of density n is 

Hi 

negligible for the range of b considered. For example, when 

—8 

a = 12 and D = 10 , an increase in b from 5 to 9 results in 

an increase in v from 320 km/sec to 340 km/sec, an increase in 

Hi 

Tpg from 2.8 x 10^ °K to 5.0 x 10^ °K while the density n £ re- 


mains constant at 11 cm 


-3 


However, since n^, remains essentially 


constant while v TC and Tp E increase (as b increases) • along ! a curve 


1/2 


of constant a, we note that at a given point (T „ , v_) in 

pHi Hi 

Figure 1 there will be a slight increase in the corresponding 


value of n_,. 

Hi 


In many treatments of the solar wind [ Parker , 1963, 1965a, 
1969; Noble and Scarf, 1963; Brandt et al. , 1965] it is usually 


concluded that a heat source is necessary between 1R 0 and 2R 0 
in order that derived coronal density distributions in this 
region can be made to fit the observed distributions. In this 

0 O 

case, the base temperature T ^ = 1.2 x 10 K at r = 2R . used 

po o o’ 

in the above, may be too low. The truth of this, of course, is 
dependent upon the base temperature at 1R 0 and the degree of 
heating required between 1R 0 and 2R 0 . We have examined the ef- 
fects of increasing the base temperature Tp Q at 2R 0 fQ? values 
of a, b, and D q corresponding to solution points of Figure 1 in 
the range of observation. In each case, an increase in Tp Q from 

1.2 x 10® °K to 1.5 x 10® °K (the value of T ) resulted in an 

eo 

increase in the density n while the remaining parameters v r , , 

Hi E 

Tp E> and T eE remained essentially unchanged. For example, 
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increasing Tp Q from 1 . 2 x 10^ °K to 1 . 5 x 10^ °K with a = 12, 

—8 

b » 5, and D = 1.5 x 10 leads to an increase in n from 

O £j 

—3 —3 

11 cm to 16 cm while v„ and T _ remain constant at 340 km/sec 

i!i P-Cj 

4 o 

and 4.5 x 10 K, respectively, Therefore, the essential change 

6 o 

required in Figure 1 when the base temperature is 1.5 x 10 K 

6 o 

instead of 1.2 x 10 K is an increase in the indicated values 
of the density. 

We have thus shown that the conclusion above stating that a 
primary heat source must be placed within about 25R 0 to attain 
wind speeds greater than 300 km/sec remains unchanged when b and 
T are allowed to vary over the above mentioned ranges. The in- 

Jr 

1 /2 

crease in density n^ at a given point (Tp E , v^,) , resulting from 

an increase in b or T can be reduced by decreasing the base 

density n Q . We find that a decrease in n Q , for a given a and 

D , results in an increase in v and T „ with a corresponding 
o hj pJcj 

decrease in n E in the range of observations. As occured when b 

1/2 

was varied, the solution points (Tp E , v E ) approximately follow 

curves of constant a. The increase in T pE as n Q is reduced is 

contrary to the result obtained in the model without heating 

(T „ decreased) . This increase in T when a heat source is 
PE pE 

applied, is due to the fact that as n Q is decreased (and con- 
comittently decreasing n(r) , for r > 2R 0 ) there is effectively 
more energy deposited per proton, part of which appears as an 
increase in the proton temperature. 

Having established the importance of heating inside 25R Q , 
consider the second class of models displayed in Figure 2 , where 
energy deposition takes place over a finite range extending from 
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the base. We have used the same representation as in Figure 1, 

in which case the solid lines are the theoretical results and 

the dashed lines are the empirical results. As before, the 

dotted lines are the analytical results lying in the region where 

pressure isotropy is not expected to be maintained by Coulomb 

collisions. In this series, solution points move from the basic 

1 /2 

solution without heating to higher values of T pJ , and v^, along 

curves of constant b as the heating rate is increased by increasing 

D q from zero. The values of the amplitude factor D Q for the 

solutions indicated by the density values and listed in an order 

7 

starting with the HS solution are: b - 2, D Q x 10 = 0.3} 1.3; 

for b = 10, D Q x 10 8 = 1, 1.5, 3; for b = 14, D Q x 10 8 = 0.4, 

1.5, 2, 2.5; for b = 18, D Q x 10 8 = 0.2, 0.6, 1, 1.5; for b = 22 , 

D Q x 10 9 = 1, 4, 6, 8, 10; for b = 26, D Q x 10 9 = 1, 3, 4, 6. 

The essential conclusions drawn from Figure 1 can be obtained 
from Figure 2 except in the instance where heating takes place 
close to the base as indicated by the solution curve for b = 2. 

In this case, when the extent of the region of the energy deposition 
is effectively restricted to within 2R Q from the base, we note that 
not only are the proton temperatures much lower than observed but 
that there is an undesirable increase in the density above that 
of the basic model without heating. This density increase occurs 
because the heating is not sufficiently extended to yield the 
higher wind speeds associated with lower densities. 

Then, increasing the extent of the heat source from the base, 
depleted by the remaining solutions in Figure 2, results in a 
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decrease in n and an increase in the corresponding values of 
h 

v_ and T _ when D is increased. As with the first class of 
E pE o 

models, there is a simple relation between the parameters of the 
heat source for corresponding solution points falling on data 


points; namely, D Q varies as the inverse fourth power of b over 
the range 275 km/sec < v £ < 400 km/sec. Densities n £ corresponding 


to V E and T pE 


in the range of observation are somewhat higher in 


the second class of models, pointing to the richness in the number 


of possible heating models. In addition, variations in T _, 

pilj 

T v„ , and n„ resulting from variations in n^ and T „ are 

ehi Jcj ill o po 

qualitatively similar to those occuring in the first class of 


models . 


To further elucidate the preceding results, let us consider 
a few examples of temperature, density, and speed profiles which 
are consistent with the empirical results. Due to obvious limit- 
ations, we present only a few typical cases of the second class 
of models (heating over a finite range extending from the base). 

The conclusions drawn from this class will be applicable to the 

first class. Specif icially, we consider the two cases in Figure 

—9 4 o 

2 corresponding to: a=2, b=26, D Q = 3 x 10 , T pE = 2.7 x 10 K, 

T g g = 3.4 x 10^ °K, v £ = 300 km/sec, n^ = 13 cm~^; a = 2, b = 22, 

D q = 6 x 10“ 9 , T pE = 3.7 x 10 4 °K, T gE = 3.4 x 10 5 °K , v £ - 320 
km/sec, n^ = 13 cm 

The electron and proton temperatures for these models, as 
functions of heleocentric radii, normalized to one solar radius, 
are shown by solid lines in Figure 3, For comparison, the results 
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of the principle model of HS without heating are indicated by 

dashed lines (as will be done hereafter) . The new results for 

the electron temperature T g have not been shown since they are 

only slightly lower than found in the principle model (e.g., 

5 o 

T gE = 3.4 x 10 K in both cases); the lower temperatures are 
due to an increase in expansive cooling above that of the basic 
model. The proton temperatures in both models are elevated 
above those of the basic model, in direct response to the heat 
sources extending from the base. In each case, the proton 
temperature declines monotonically , although becoming almost 
isothermal around IORq. (In the models of the first class, 
the temperature profiles are not always monotonic due to their 
response to the position and shape of the heat sources •,) Of 
the two examples shown, that corresponding to the upper proton 
temperature curve lead to the higher wind speed at 1 AU due to 
the larger overall pressure gradient (the corresponding densities 
are about the same) , resulting in greater acceleration. These 
temperature profiles serve to remind us that the higher tempera- 
tures T „ and wind speeds v,,, agreeing with observation, correspond 
pE E 

to overall energy deposition closer to the sun; i.e., as Tp E 

and v increase, the half width b decreases (and correspondingly D 

increases) . 

In the basic model without heating, proton thermal conduction 
was essentially unimportant beyond a few solar radii from the 
base. However, when heating takes place, proton thermal conduction 
is important to greater distances from the base due to the elevated 
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proton temperature along with the strong temperature dependence 
of the conduction coefficient. In particular, for the examples 
shown in Figure 3, proton thermal conduction accounted for about 
20% of the proton energy flux out to 25-3OR 0 from the base. 

Beyond this point it becomes unimportant and only thermal 
convection and expansive cooling are important. At several 
hundred astronomical units and beyond, thermal convection, ex- 
pansion cooling, and energy exchange are dominant processes and 

/2 /n 

T “* r ' as r - 
P 

Figure 4 gives the density and speed profiles for the two 
examples shown in Figure 3. As noted previously, the example 
with the highest wind speed corresponds to the case of highest 
energy deposition leading to the greatest pressure gradiant 
and acceleration near the base. The density distributions for 
both cases are about the same, although somewhat lower than that 
of the basic model. It is worth noting in Figures 1 and 2 that 
the densities n„, lying in the range of observation with 

ili 

Vg ~ 300 km/sec, are essentially constant for a given class of 
models. The density decrease near the base, relative to the 
basic model, is permissible as can be noted in Figure 5 where 
the theoretical density distribution is compared with coronal 
density observations at solar minimum reported by Blackwell 1*19563 , 
Michard [1954] , and Allen [1963] and the representative dis- 
tribution by Newkirk [1967] . 

The effects of raising the base proton temperature Tp Q 
are shown in Figure 6 and 7, where the case corresponding to 
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D Q = 3 x 10 , a = 2, and b = 26 is considered. The dot- 

6 o 

dashed lines in these figures result when T = 1.5 x 10 K 

po 

while, for comparison, the solid lines are based on 
6 o 

T = 1.2 x 10 K. The increase in proton temperature at the 

ir 

base, resulting in an increase in the pressure gradient near 
the base, leads to a higher base wind speed and corresponding 
larger particle flux. However, as mentioned above, the proton 
and electron temperatures and wind speed are essentially un- 

Q 

changed at 1 AU while the density is increased to 20 cm (to 

maintain constant particle flux) . The model with higher particle 

flux has lower acceleration near the base since more particles 

must be accelerated with less energy available per particle. 

The undesirable high density at 1 AU can be reduced by decreasing 

5 -3 

the base density. For example, a base density n Q = 8.5 x 10 cm 

leads to n E = 10 cm -3 , with v £ = 330 km/sec, T pE = 4.9 x 10^ °K, 

5 o 

and T gE = 3.4 x 10 K. Much further density decreases at the 
base leads to poor fits to the observed coronal density distri- 
butions (2-20R Q ) . 

In the preceeding, the effects of viscosity were ignored. 

Due to low proton temperatures, viscosity is negligible in the 

HS model without heating. However, when heating takes place, 

viscous terms become larger as the proton temperature rises -since 

5 /2 

the coefficient of viscosity varies as T .We have estimated 

P 

the effects of adding viscous terms [see Braginskii , 1965] 
to the momentum and heat equations by treating them as small 
corrections and deriving their values from corresponding preceeding 
temperature and velocity profiles [n-th profiles; see Hartle 
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and Sturrock, 1968] in the iterative sequence. We find that 


collisional viscosity leads to small corrections for those 
cases corresponding to the upper portions of the solid lines 
(v > \> eX p) in Figures 1 and 2. For example, considering 
the upper ends of the solid lines, we find that for the range 
Vg ~ 275 km/sec to ~ 375 km/sec the proton temperature in- 
creases ~10% to ~15% and v E decreases slightly (£ 2%) when 
the viscous terms are included (a consequence of the fact 
that viscous heating is greatest for r £ 30 R Q ) . Below the 
endpoints of the solid lines, towards lower temperatures, the 
viscous corrections diminish. For those solution points with 
Vp < v eX p (dotted lines) in the observational range, correspond- 
ing viscous corrections are somewhat larger and result in an 
increase in Tp E by a factor of about 2-3 when v E is high, 

350-375 km/sec. However, Parker [1965a] has noted that when 

collisions are rare the viscous coefficient should be reduced 

2 

by a factor (r/\) , in terms of the collision mean-free path 
hence, the above mentioned increase in Tp E by a factor of about 
2-3 is an overestimate. With these considerations, we note 
that viscous corrections are not large over the parameter range 
of interest and that conclusions made above are essentially 
unaltered . 
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4. SUMMARY AND CONCLUSIONS 


We have argued, as others have, that observed solar wind 
speeds and proton temperatures above those of HS suggest the 
need for additional energy transport across the base of the model, 
possibly in the form of waves which then deposit their energy in 
the system by means of collisionless damping. To account for 
such nonthermal energy deposition, we have taken a phenomeno- 
logical approach by adding a<i hoc a heat source to the proton 
heat equation without relating it to any particular mechanism. An 
electron heat source was not included. For ease of analysis and 
interpretation, the form of the proton heat source used was a 
density weighted gaussian given by equation (11) l ' The batsas©f'the 
model was taken at 2R 0 where corresponding values of temperature 
and density obtained by H&. were used as boundary conditions. 

The empirical relation between v„ and T „ determined by Burlaga 

E pE — 

and Ogilvie Cl970] was used as a constraint on the parameters 
describing the heat source. 

In every case considered, we found that, when the amplitude 

of the heat source with a given width and position was increased 

from zero, the wind speed v and proton temperature T „ increased 

E pE 

while the electron temperature T „ decreased slightly from their 

respective values obtained in the HS model without heating (the 

estimates of v £ and Tp E in the HS model are in good accord with 

the corresponding values determined from the empirical relation 

of Burlaga and Ogilvie). In most instances, when the intensity 

of the heat source was raised, the density n„ was found to decrease 

E 
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_3 

from 15 cm obtained in the HS model; however, when energy 
deposition took place essentially within 2R 0 from the base, 

_3 

the density increased from 15 cm when the intensity of the 
heat source increased. 

More specifically, upon comparing the variations of v_ and 

ill 

T „ with the corresponding results from the empirical relation 

pHi 

of Burlaga and Ogilvie, we find that an increase in the energy 
deposition rate, with deposition taking place beyond about 30 R q 
from the base, results in a significant increase in T ^ from 
the HS result and a negligible increase in v . In this case, 

£j 

upon increasing the heating rate, n does not decrease more than 

hi 

_3 

7% from the 15 cm of HS. In the opposite extreme, relative 

to the results of Burlaga and Ogilvie, v^, increases substantially 

while T „ rises slightly when the heating rate is increased in 
P-Ei 

sources depositing energy within 2R Q from the base. With such 
an increase in the heating rate the corresponding density n„ 

hi 

-3 

increases more than 30% above 15 cm . However, when the primary 

energy deposition takes place within the range of approximately 

2R^ to 25R_ from the base we showed that v_ and T _ can be brought 

into direct correspondence with the empirical results of Burlaga 

and Ogilvie (over the range 275 km/sec < v„ < 400 km/sec, 10^ °K < T 

k pE 

4 o 

<8x10 K) when suitable variations in the parameters of the heat 

source are made. In particular, refering to the heat source of 

equation . (11), such correspondence obtains when D Q « b~^ with 

—2 

a = 2, and D Q « a with b = 5 ( of course, other parametric 
relations are possible) . The corresponding proton density at 
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-3 -3 

Earth’s orbit ranged from 11 cm to 13 cm , a reduction of 
12% to 25% below the value obtained by HS. Such a reduction is 

_3 

desirable since the observed average density is about 5-6 cm ; 

_3 

however, the observed density is known to range from 1 cm to 

o 

30 cm [ Neugebauer and Snyder , 1966b; Wolfe et al , 1966; 

Lazarus et al , 1966; Gringauz , 1966]. 

Based on these considerations, we concluded, in the context 
of a two-fluid, spherically symmetric solar wind model, that 
primary energy deposition should take place inside about 25R Q 
from the base to obtain the more commonly observed wind speeds 
v E ~ 275 km/sec while simultaneously maintaining corresponding 
values of T p low enough to agree with data. Parker [1963, 1969] 

P-Ci 

and Burlaga and Ogilvie [l970] have also suggested the need for 

an energy source in this region. This result is consistent with 

the nonthermal heating models recently proposed by Barnes 

[1969], and Fredricks [1969]. On the one hand, Barnes has proposed 

that the collisionless damping of fast-mode hydromagnetic waves 

propagating outward from below the base will deposit most of their 

energy into the proton (and electron) gas within 25R Q . On the 

other hand, the turbulent ion-acoustic wave heating model proposed 

by Fredricks is operative over the range 20R G to 60 Rq when the 

conditions of the HS model prevail. In the Fredricks model ion- 

acoustic waves are generated (becoming unstable when T /T & 7-10) 

e p 

by compressive MHD perturbations which could originate in part 
below the base (thereby constituting an external source) . Based 
on the preceeding, the proton heating mechanism proposed by Barnes 
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is expected to be more significant since it >is operative over more 

of the important range 2R Q ~ r ~ 25R Q while that of Fredricks 

is only operative over the outer end of this region. In addition, 

after some heating has occurred, the inner boundary of the Fredricks 

model, r ~ 20R^, will move out further to satisfy T /T ~ 7-10 
© e p 

while proton heating in the Barnes model over the range 2R 0 ~ 

< 2 
r ~ 25R 0 is expected to be more effective when 8 = 8tt n k T /B 

is increased. 

Even though we have shown that it is possible to obtain ac- 
cord with the empirical relation between v„ and T _ by application 

£< Pill 

of energy sources over the region 2R Q ~ r ~ 25R 0 , we wish to make 

it clear that we do not rule out the possible importance of heating 

over the range r ^ 25R 0 . For example, consider the solid curves 

in Figures 1 and 2 lying well below the empirical results of 

Burlaga and Ogilvie. It can be shown that solution points on 

these curves can be raised essentially vertically to obtain 

correspondence with observation by an additional heat source 

depositing energy anywhere over the range 0.12 .AU ~ r ^ 1 AU. 

Of course, the bulk of the heating still takes place inside 25R 0 . 

For such an "extended” or "double -humped” heating model the 

corresponding densities n E do not change noticably from their 

original, somewhat high values as the solution points are raised 

upward towards 'higher Tp^when a; heat source -is applied beyond 25R 0 . 

If a small amount of heating beyond 25R 0 (^1-5% of total 

energy) is actually required to raise T _ , then the turbulent 

P-Ci 

ion-acoustic heating model of Fredricks, also operative in this 
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region, may be important. In addition, two other possible external 
heating mechanisms in this region have been proposed by Coleman 
[1968] and Jokipii and Davis [1969]. Jokipii and Davis pointed 
out that the observed hydromagnetic wavelengths ~ 10 km are 
probably generated in the solar wind by the "collision" of fast 
and slow streams. The hydromagnetic waves generated then heat 
the protons locally through the collisionless damping mechanism 
of Barnes [1966] . We point out that such a hydromagnetic heating 
process in the very low P + streams of the HS model would be re- 
latively insufficient and therefore the cooler streams should 
arrive appropriately pre-heated at the point of "collision" for 
the Jokipii and Davis model to be of significance. Coleman, on 
the other hand, has suggested that energy extracted from differential 
streaming motion cascades down through a hierarchy of turbulent 
eddies to those approximately the dimensions of a proton Larmor 
radius at which point dissipation takes place by proton cyclotron 
heating. However, Burlaga and Ogilvie [1970] have observed local 
heating or "hot spots" in the vicinity of positive velocity gradients 
but not negative gradients and therefore conclude that the heating is 
not a result of velocity shear but due to the collision of fast 
and slow streams. They also noted that the velocity gradients 
and corresponding "hot spots" occur infrequently and therefore 
the associated heating mechanism is not a dominant large scale 
phenomenon. 

Proton heating in the region r ~ 25R Q by means of internally 
generated unstable waves has been proposed by several investigators 
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[ Parker , 1958 a, b, 1963; Scarf et al. , 1967; Kennel and Scarf , 
1968; Hundhausen , 1968; Forslund , 1970]. Since only a small 
fraction of the total energy is required to raise the proton 
temperature in this region, solution points in Figures 1 and 2 
lying below the empirical results can also be raised approximately 
vertically by the application of an internally generated heat 
source; i.e., the corresponding v„, n and T _ will not change 

i !i ili pili 

appreciably from the case when an equivalent external heat source 

is applied. It is clear that we can not point to any factor or 

factors in the present analysis that would rule out the need for 

heating beyond 25 Rq or, if such heating is necessary, whether 

or not it is derived from an internally or externally generated 

source or both. We only argue that it is possible to account 

for the observed relation between v_, and T by application of an 

Hi PUj 

external energy source inside the heleocentric radius of about 
25R 0 and that the corresponding densities n^ are desirably lower 
than those resulting from an "extended" or "double -humped" heat 
source . 

In addition to the above internal sources operative over 
r ~ 25R q , Forslund Cl970] has shown that an internally generated 
proton heat source is possible inside 25R Q . He pointed out that 
the expected skewing of the electron distribution, due to the 
high rate of thermal conduction, can be responsible for exciting 
unstable ion-acoustic waves leading to turbulence from which 
protons can be heated through linear and non-linear Landau damping 
and collisional damping of the waves oa the ions. The inner 


32 



boundary for this mechanism occurs where T e /T p ^ 7, which can 

take place as close as r ^ 10R^ when the conditions of the HS 

model prevail. However, when heating occurs, the inner boundary 

will move outward to satisfy T /T £ 7. Forslund has also 

e p 

pointed out that such ion-acoustic turbulence should lead to a 
reduction in the thermal conduction coefficient. 

We have argued that energy must be added to the system to 
account for the observed high speed streams. We do not expect 
that a reduction in the electron thermal conduction coefficient 
at r ~ IORq, as suggested by Forslund, can effect an increase 
in the total energy flux. It seems unlikely that more energy 
will be drawn out through the base, say, by a possible increase 
j in conduction flux at the base, when the conductive coefficient 
^ is reduced at r ~ 10R Q . However, the relative merits of external 
and internal energy sources can not be conclusively evaluated 
until the macroscopic effects of the latter have been quantita- 
tively investigated. 
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Fig. 1 - The square root of proton temperature T E at 1AU 

in kilo-degrees Kelvin versus wind speei vg in 
km/sec. Analytical results are given by solid lines 
(and corresponding dotted lines) with corresponding 
values of density ng in cm-3 indicated at selected 
points. For ' comparison , the empirical results of 
Burlaga and Ogilvie are shown by dashed lines. 
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Fig. 2 - The square root of proton temperature Tpg at 1AU 
in kilo-degrees Kelvin versus wind speed vj; in 
km/sec. Analytical results are given by solid lines 
(and corresponding dotted lines) with corresponding 
values of density ng in cm -3 indicated at selected 
points . For comparison, the empirical results of 
Bur lag a and Ogilvie are shown by dashed lines . 




Fig . 3 - Proton temperature T p and electron temperature T e 
versus radius r. Temperature profiles resulting 
from heat sources characterized by a=2, b=22 , 

D 0 =6 x 10 -9 , and a=2, b=26 , D 0 =3 x 10“9 are given 
by solid lines. For comparison, corresponding re 
suits of HS model are shown by dashed lines. 
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Fig. 4 - Wind speed v and density n versus radius r. Density 
and wind speed profiles resulting from heat sources 
characterized by a=2 , b=22 , D 0 =6 x 10 -9 , and a=2 , 
b=26 , D 0 =3 x 10 -9 are given by solid lines. For 
comparison, corresponding results of HS model are 
shown by dashed lines. 
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Fig. 5 - Density n versus radius r. The solid line corresponds 
to the density profile of Fig. 3, given by a solid 
line. The observational results of Blackwell, Michard, 
Allen and Newkirk are indicated. For comparison, 
the density profiles of HS model are shown by dashed 
lines . 




Fig. 6 - Proton temperature Tp and electron temperature T e 

versus radius r. Solid and dot-dashed lines correspond 
to a heat source characterized by a=2 , b=26, and 
D 0 =3 x 10“9 . Solid line corresponds to T po =1.2 x 10® 
and dot-dashed line to T po =1.5 x 10 6 °K. ror com- 
parison, corresponding results of HS model are shown 
by dashed lines. 





Fig. 7 - Wind speed v and density n versus radius r. Solid 
and dot-dashed lines correspond to a heat source 
characterized by a=2 , b=26, and D 0 = 3 x 10 - 9 . 

Solid lines correspond to T po = 1.2 x 10 6 °K and 
dot-dashed lines to T po = 1.5 x 106 or. p or com- 
parison, corresponding results of HS model are shown 
by dashed lines. 




